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Abstract. We study an homogeneous irreducible markovian random walk in Z", with an 
antisymmetric perturbation acting only in one point. We compute exactly spatial correction to 
the diffusive behaviour in the asympotics of probability, in the spirit of local limit theorems for 
random walks. 

1. Introduction, definitions and main result 
In this work we are concerned about a Markov random walk in Z 1 ', with transition probabilities: 
(1) U(x,x + u) = P(u) + 5 xfi c{u). 

Here P is the probability of a symmetric homogeneous walk, and the other term is a perturbation 
acting only in the origin. This is equivalent to a free random walk in which only the transition 
rates from the origin are modified. We are interested in finding the probability Ii n {x) of being 
at point x at time n, when n — > oo. 

It is well known that for the homogeneous case a local limit theorem holds: the probability 
distribution P n (Q,x) of being in x starting at the origin is asymptotically of the form: 

P n (0,x)= — j - axpf- Cg^T ) [ 1 + y_LQ j (^]\ +rn{x)o f 1 



\B\(2im) v / 2 




n {v+r-2)/2 



where P n (0,x) is supposed to have the first r > 3 moments finite, B is the covariant matrix of 
the process (we have supposed no drift) and Qj are certain polinomia of degree 3j. Moreover, 
r n (x) is a correction for small x: 

\x\ 



For a more detailed discussion see [JJ. It is important to stress the range in which this classical 
result holds: of course it is meaningful only when the magnitude of the corrections is not bigger 

2 

than the leading term; in this case up to the scale x = 0(nz). 

Despite the rather simple formulation of the problem, it seems that no extensions of the previous 
theorem exist in the case in which the probability rate is changed only in one point of the lattice. 
However the obstacle problem also arises in the setting of random walk in random environment, 
see for instance [I], [9], [2]. 

In [6] is proven a general result about a random walk with an homogeneous part plus a per- 
turbation V(x, y) acting in a finite neighborhood of the origin. In this case a correction to the 
diffusive behaviour arises. If we set 

A n (0,a:) = n n (0,x) -P n (0,x), 
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we have for v = 1 

1 _ ( 1 

A„(0,x) = — p= sign(x)e z^ 1 + ip n (x) + o 



>n " ' \n v l 2 
with ip n {x) decreasing at least polinomially, and for v > 2, 

1 C ( 1 

A «^°' ^ ~ 1 + | x |(i/-l)/2 + 

This result is achieved under the general assumption that the perturbed random walk is irre- 
ducible (see [8]) and that there are no traps {i.e. a finite set TcS" such that a walk starting in 
T will never leave this set with probability one). Furthermore in dimension v = 2 an additional 
and more technical hypothesis is needed, for which we refer the reader to the original paper [6]. 
Anyway the method used there is based on sophisticated functional analysis tecniques, and so, 
as the same authors write, a direct probabilistic intepretation misses. 
Other interesting works in the subject are [7], |10] . 

In our work we analyze the problem in the simpler case in which the perturbation is concentrated 
in the origin, but by a method with a clear probabilistic interpretation, that is the Gihman- 
Skorohod expansion of the characteristic function of the random walk. 

In order to clarify our argument, we will separate the symmetric and antisymmetric part of the 
potential function, c(u) = es(u) + a(u), for a certain parameter e > ruling the strenght of the 
symmetric part of the interaction. 

Moreover we assume P has only the first three moments finite, and c the first two. Anyway, in 
principle we could relax this assumption, allowing more regularity both for P and c; but it will 
be clear in a moment that for our purpose it is not useful at all, since at the diffusive scale higher 
moments play no role. Actually, since we will deal with P a symmetric and c an antisymmetric 
function, we will need only the covariance of the free walk and the drift given by the perturbation. 
Finally, we define a brownian scale, x ~ y/n, when x = 0(y/ncp n ), with (ft n an increasing 
sequence with n growing slower than any power. For example, the result in [6] is valid up to 
x = 0(\/n logn). 

In Section 2 we find the reprensentation formula for the transition probabilities. It turns out 
that relevant simplifications occur by setting e = 0, so in Section 3 we focus only on the anti- 
symmetric case, and we find the asymptotics. Performing the same calculation in the symmetric 
case is harder from a technical point of view, due to the more involved representation formula. 
Neverthless we believe that the method we use can be applied straightforwardly to this case too. 
Our main result is the following: 

Theorem. Be given an irreducible markovian random walk in TLf , with transition probabilities 
given by (CP, with the perturbation antisymmetric, and 

P{x) = 1; xkp ( x ) = °> fc = 1, 3; ^2 x 2 P(x) = B < oo; 

and 

xc[x) = d < oo; ^2 x k c{x) = 0, k = 0, 2. 

Then the following asymptotic form for the probability of being in the point x at time n up to the 
scale x ~ y/n holds: 

n »(°' I » = oW^ exp (- I ^) +A " (l) ' 
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where we have the following form for A^(x); 



Kix) = \d\ cos(d, x) E 5 v n {x) +o. r , ., 

rj2 \n u / z 



with 



1 



l- 



i 



(2) 5n{x) ~ \B\^ n»l* Ji da a^(l - aW 2 

n 

is a 0(1) function of n when x is fixed. 

Remark 1. We will prove the following properties for 6 n (x): 

(1) in dimension v > 1, 5^(x) is O(l) in n when x is 0(1) in n, i.e. near the perturbation; 

(2) in dimension v = \, S^(x) gives a relavant contribution also for x growing at the browian 
scale x ~ \fn. 



Remark 2. Although we use the same tecnique exposed m|, we get that our theorem holds for 
a smaller range of values of x, i.e. up to the scale \fn (essentially as in jBJJ rather than ra 2 / 3 . 
We believe that this is not a crucial feature of the problem, but it is rather due to the fact that we 
focus on the main term, and we estimate the higher order correction roughly with respect to the 

more careful analysis made in [3]. Namely we let the term o ^^72^ of unspecified form. For the 

same reason, for example, we have that the regularity of the probability, i.e. the number of finite 
moments of P and c, play substantially no role in our proof, while of course it will be relevant in 

the exact form of the remainder o ( — ^j. 



Recently, a very similar theorem has been obtained with different methods for v = 1, but with 
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a perturbation having no definite parity [3], up to the scale x ~ n?. We stress that the two 
approaches bring to analogous results. 

2. The representation formula for the transition probability 

The aim of this section is finding a representation for the transition probability (|TJ) . 
It turns out that it is useful to define the probabilities II„(a;, y) and Pn{x, y) of going to y starting 
from x, without passing through the origin, respectively of the perturbed and free walk. We will 
concern expecially about P , and we can represent such probability in terms of P, by using a 
classical argument, namely the inclusion-exclusion principle, as stated by the following 

Proposition 1. We have 

(3) P n °(0,x) =P n (0,x)-p n (x), 
where p n (x) is a positive function of x for every n: 

n-l 

(4) Pn ( x ) = Y J {-l) p Pn-kM0)P kl ^ k2 (0,0)...P kp _ 1 - kp (0,0)P kp (0,x) 

p=l fcl>fc 2 >...>fcp>l 

Proof. The unperturbed probability to reach a point x from the origin can be decomposed by 
the following expression: 

n-l 

P n (0,x) = ^P n _ fc (0,0)P fe °(0,x), 

k=l 
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i.e. a sum on all the possible cicles around the origin, finally followed by a direct path to x. 
It is straightforward to invert the previous formula in order to recover (|3|), just noticing that 
P kx (0,x) = P% x (0,x), where 1 < k x = mm{k : P k (0,x) > 0}. □ 

Now we focus only on the antisymmetric part of the perturbation to obtain in primis a formula 
for returns in zero, as can be seen by the following 

Proposition 2. Let IT stated and P be defined as in f7]j, but considering only the antisymmetric 
part of the point potential i.e. e = and c(u) = a(u). The following representation formula 
holds: 

n-l 

(5) n n (0,x) = P n (0,x) + ^P fc (0,0)(a*P n °_ fc _ 1 (x)). 

k=0 

Proof. We have 

n°(o,o) = £n(o,i/)id(y,o) 

2/e£ 

= ^T(P(0,y) + a(y))Pti(y,0) 

= Yl p (°< y) p n-i (v, o) + E <y) p n-i (o, -v) 

= ^(0,0), 

because of the antisymmetry of a(y). Hence the probability of returning in zero for the first time 
in k steps is equal to the unperturbed one. Furthermore, for an arbitary time n, we have 

P n (0,0) = ^P fc °(0,0)P n _ fc (0,0), 

k 

n n (o,o) = ^n°(o,o)n„_ fe (o,o) = ^p fc °(o,o)n n _ fc (o,o). 

k k 

This is a recursive formula for P n (0, 0) and II n (0,0), with the same initial condition P(0, 0) = 
11(0,0), that implies P n (0,0) = II n (0, 0). Now we are ready to find a form for the probability 
ll n (0, x). We use a decomposition with respect to the time of first return at the origin, such that 

n—l n—1 
U n (0,x) = ^n fc (0,0)n°„ fc (0,x) = ^P fc (0,0) £ (P(0,y) +a(y))P°_ k _ 1 (y,x) 

k=0 k=0 y& v o 

n-l 

= ^P fe (0,0)(P n _ fc (0,^)+ Y, a(y)P^ l _ k _ 1 (0,y-x)) 
k=o yeZ"o 
n—l 

= P n (0,x) + ^P fc (0,0)(a^P n °_ fc _ 1 (x)). 

k=0 

□ 

Remark 3. It is easily seen that the formula |5|) can be simplified: infact the paths returning 
into the origin, due to the antisymmetry of a{x), give a vanishing contribution to the convolution, 
i.e. a * P®_ k _ 1 (x) = a* P n -k-i( x ) an d thus we can rewrite the previous representation formula 
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as 



n-l 



(6) n n (0,x) = P n (0,x) + ^P fc (0,0)(a*P„_ (fc „ 1 (x)). 

k=0 

Now we can concern about the symmetric part of the perturbation. In this case, due to the 
symmetric nature of the interaction, formulae easily become very complicated; for this reason in 
this work we will treat in full detail only the case e = 0. 
We can state and prove the following 

Proposition 3. Be U and P defined as before in (CP, but considering only the symmetric part 
of the one point potential, c(x) = es(x). Thus the following representation formula holds: 

n — l n— 1 I 

n n (0, x) = P n (0, x) + e ^ Pl(0, 0) (s * (a;) + e £ ]T fl-fc(0, 0)(a * P^_ 1 )(0)P n °_ / (0, x) 

1=0 1=0 k=l 

n—l u I a 

+ E ^(a)(o,o)p n %(o,x)n(^nVi)( ) 

1=0 a=2 K(a)=l 1=1 

n—l u ( a 

( ? ) + E ^(a)(o,o)( s *p n v 1 )(^n(^ p fc Vi)(°)- 

Z=0 a=l K(a)=l i=l 

Proof. Following the same scheme adopted in the proof of the previous lemma, we start by 
computing the probability for returns in the origin in n steps: 

n n (o,o)= n° fel (o,o)...n° feM (o,o) 

fci + ...+fe M =ra 

with of course fi = n/2. Again we have 

(8) II°(0,0) = P fc °(0,0) +e(s*P fe °_ 1 )(0), 
and 

(9) n°(0, x) = P fe °(0, x) + e(s * Ptx){x). 

Therefore, putting ([8]) in the representation for II n (0, 0), and naming K(a) = Y2t=i we nave 

n n (o,o) = Yl (n° 1 (o,o) + 6( s *p fc 1 _i)(o))...«(o,o) + e(^p fc ° M _ 1 )(o)) 

fci + ...+fc M =ra 

[in a 

= P n (0,Q) + Ye a Yl p n-K(a)(0,0)l[(s*I%-i)(0)- 

a=l K(a)=l i=l 
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Now we are ready to write down the representation formula for the transition probability from 
the origin to a given point: 

n—l 

n„(o,x) = ^n,(o,o)n°_,(o,x) 

1=0 

n—l n—l [i I a 

= ^PKo,o)n°_ i (o,*) + ^ J>« ^-W ' ) n ^( ' )n( s * p *Vi)(°) 

1=0 1=0 a=l K(a)=l i=l 

n—l n—l I 

= P n (0,x) + £ ^PK0,0)(5*P^_J)(x)+ £ ^^^_ jfc (0,0)( S *P fe °_ 1 )(0)P n %(0,x) 

1=0 1=0 k=l 

1=0 a=2 K(a)=l i=l 

n—l [i I a 

+ EE £l+Q E Pi-Kiamo^s+p^w^s+p^m. 

1=0 a=l K(a)=l i=l 



Putting together these three propositions, we obtain the subsequent 

Lemma 1. The following representation formula for the transition probabilities holds 

n-l 

n n (0,x) = P„(0,x) + ^P fc (0,0)(c*P n _ fc _ 1 )(x) 

fc=0 

n—l n—k 

+ e^^P„-/ l - A; (0,0)( S *P h _ 1 )(0)P fc (0,x) 

k=0 h=l 

n—l u I a 

+ EE £ ° E fl-jf(«)(0,0)n( a *i^_ 1 )(0)J*_ I (0 > x) 

/=0 a=2 K(a)=l i=l 
n—l u I a 

+ EE £l+Q E ^( Q )(o,o)n(^nVi)( )( s ^nVi)w 

/=0 a=l K(a)=l i=l 

(10) - R n (x), 
where 

n-l i 

R n (x) = e^^P i _ fc (0,0)( S *P fc °„ 1 )(0)p n _ / (0,x) 

1=0 k=l 

n—l u I a 

+ EE e ° E P i~K ia )(0,0)H(s*P^ 1 )(0)p n „ l (0,x) 

1=0 a=2 K(a)=l i=l 
n—l [i I a 

(11) + EE £l+Q E ^(«)(o > o)n( fl *i^_ 1 )(o)( fl *^_,_i)(x). 

1=0 a=l K(a)=l i=l 



□ 
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3. Proof of the Theorem 

Hereafter we will set e = 0, so that the representation formula for the transition probability 
reduces to ([6]). 

In order to prove the Theorem, it is convenient to split the proof in two steps. The first be 
essentially a convenient way to write the inhomogeneous term of the transition probability, 
following the method in [4]. At first it is useful to define the return probability in the perturbed 
point at time step n in dimension v as 

p" n = P n (Q,0). 

Remark 4. Throughout the whole paper we will use the well known behaviour of the returns 
probabilities of the free walk p v n ~ n~ v l 2 [8] . 

Therefore we are ready to state the sequent crucial lemma: 

Lemma 2. The correction term A n (x) at the leading order in n has the form 

n— 1 v , 

* — "\ V l i (B~ L x,x) 

(12) A n (x) = £ P -^^ e --^{d ■ x), 

k=i k 2 

where d = J2 X xa ( x )- 

Remark 5. Before the proof, it must be noticed that the order in time of A n is the right one, 

_v_ (B~ 1 x,x) 

that it n 2 . This is easily seen for instance by an expansion of the function e 2k 



) _ (ff-^x) 1 {B- X x,xf _ 1 (B^ x.r 



3 



2k + 2 Ak 2 6 8k 3 + 



and, for v > 0, we have that for every q > [5] 

sr^ Pn-k-i _ Const 

v + 2 + q - u / 2 ■ 

k=l k 2 11 

Proof. In primis we calculate the characteristic function of the transition probabilities starting 
by a rearrangement of formula 

n-1 

n (A) = 0°(A) + ^K- fc -i«(A)^(A), 

k=l 

where A G [— 7T, 7r] y , </> n (A) is the characteristic function of the free random walk, and 

2(A) = < x Y Xx 

P k = J2 Pk(0,x)e iXx . 

Then we write 

1^(0, x) = / #47^Vn(A) = / ^e-^ x MX), 

and we will focus on the inhomogeneous part of the characteristic function, and on the sum of 
integrals 

n— 1 

(13) A n (0,x) = Y.Pn-k-i [ ^-e-^a(\)P k {\). 
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The main idea is very simple [4j: to perform a Taylor expansion of the Fourier transform of the 
potential, that is supposed to be C 2 , and of the free walk transition probability, that has three 
moments finite. It is 

(14) a(\)=i(d-\)+r a (\) 

where d is a kind of drift term given by the first moment of a(x) and 

P(X) k = (l-I (j BA,A) + rp(A)) fc 

= e fclog(l-|(BA,A)+rp(A)) 
= e -|(BA,A)-fer P (A)) ) 

where B is the covariance matrix of P. Moreover the reminders r a (A),rp(A) are bounded func- 
tions such that 

hm = 0. 

A^O |A| 3 



We split the quantity of interest (|13p into two pieces: 



/ = / ^-e-^e-^\{d-X)e~^ 



fc=l 
n-1 



11 = E^-fc-i / 



d" A 

(27T)"' 



-ia:A e -f(BA,A) 



o(|A| 2 )e-^ (A) . 



We define also J and II as /, II with the integral taken in the complement to W of the cube. 
Now we have to notice two crucial features: 

(1) Inside [— ir, tt] u , e~ krp ^ can be replaced by 1 with an error o(l); 

(2) I and II are exponentially vanishing. 

To prove point 1 is useful to rescale A — > \y/n in such a way the integration domain becomes a 
cube of size A= and observe that 

rp(A) ~ f -7= ) */>(A,n), ip(\, n) — > if n — > 00, 



so 

e -fcnf rp(A) „ g -^(A,n) = j + q(i) [n ^ 



>n 

at most, i.e. when k is finite, while if k ~ n in the sum the terms are exponentially vanishing 
because of the presence of e~2^ BX ' X \ 

The proof of point 2 is a bit longer. It is useful to rescale the modes: A — > ^= in such a way the 

integration domain becomes a cube of size yjn, [—i/nir, yjnn} v . We will focus on the lattice of 
translations: 

[— \/wr, \/nir] u + Aq, Aq = h2ir\/n, h G 
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This of course is a partition of W, so 



/(•)=/ (•)+£/ 



<[-^/nn,^/nn] v Xq ^. q J [- ^/nn, y/nir]" +Xo 

Now we expand a and P around each Ao- It is 

a(X) = a(A ) + id(\ ) ■ A + r a (A), 



with 



a(A ) = J>* Aoa; a(*), 
d(A ) = Yl e iXt>x xa(x), 



and 
where 



P fc (A) = e klogP(\ ) e -%(B(\ )\,\) e -kr P (\) 

P(A ) = J> lAo ^P(x), 



B(Ao) = 



^(Ao) 

1 = EK-.-i £ / 7^;e^ e -^^^(a(Ao) + ^(A ) • A) 



So we can write for instance for /: 

rt— 1 



k=l Ao^O 

and it is immediate to show that 

n-l 



7 - E^-iE e " 



logP(Ao) 



fc=l Ao^O 
n-l 



E^-,-1 E e-f^^)e-^(Ao) 

fc=l A ^0 

< EK- fe -i E e - (BAo < Ao) . 

fc=l |A |>27r A /n 

where the main contribution in the sum is given by small k, and we have the same superexpo- 
nential decay as in a gaussian tail starting by 2iTy/n. 
Of course the same is valid for II. 

Moreover, we claim that even II gives a correction of small order with respect to the diffusive 
scale. Infact inside the rescaled cube we have 

Vn 

n— 1 

fc=i J 7s ^ > 

On the other hand we have seen that we can extend the integration to the whole space with a 
very little error. So we have 
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II 

T 



Z^fc=l Pn-k-1 Jr^ e e 



iiA -|(BA,A)| \ 1 2 



|Apo(l) 



>A»ELiiC*-i/R» e " Ae 



ixXp— 2 (-B^j^) 



(d-A) 



1 \ Efc=lPn-fe-l e 



(B~ x :e,:e) 

2fc M 



(15) 



\x o 



So, since it is easily seen that I = O ( — ^ ) , we have that II = o [ } +1 I if we choose x up to 



the scale 0(y/ri(j) n ), with n a sequence increasing with n slower than any power. 
Therefore we say that the mean contribution is given only by: 

n-l 



(16) 

That leads to 



fc=l 



A n (x) 



fc-i 



d"A 



u (2tt) u 

n-1 

Gi"' 2 ^ 
1 fc=i 



ix\ -\ (BA,A) 



i(rf-A). 



™ _ i v i 

1 X ^ Pn-fc-1 -^-'■■■) , , v 

v+2 e 2k (d-x). 
k 2 



□ 



Now we can rewrite the correction (fT2l) as 



(17) A n (x) = l*^^g%l e -^. 

I«l 7 fc=i fc 2 

Since we are interested in the asymptotic behaviour of the correction term, we can replace k _-y 
with (n — k — X) v l 2 . With further manipulations we get 

n— 1 (B~ 1 x,x) 

\x\ \d\ cos(d, x) 1 ^ e 2 ™ a 
= FbKV^ n^a^/2+i(l-a)^/ 2 

_ |x||rf|cos(ri,z) 1 e~ 2 "°' fj_ 

n 

where we have set a = fc/n, and we have written the dense sum in [0, 1] as an appropriate integral 
plus a small error. So we have 

(b~ 1 x,x) . _ \ 

(19) A n (x) = |d|cos(d,a) 5 n (x) + o (— j , 

with 

I I -, l-c (B~ 1 x,x) 

\X\ I n e a 2n 

(20) 5 n (x) = -Ut^-t, / da 



B \v/2 n v/2j L a v/2+U l _ a y/2- 

71 V / 

This concludes the proof of the theorem. 
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For sake of completeness we analyse the integral in (fl8|) separately in the case of dimension 
v = 1, 2, 3 in order to get exact formulas for the asymptotics. 

Case v = 1: 

In one dimension we put B = a 2 and the integral turns out to be [5j 

1 ?. 2 / / / x ' x x 

\/27rn * 2 / / i n — 1 



e v Z7rn * / r / F / n — 1 I r / Fl 

da — : = = — : — : — e 2n CT 2 u e rf _ a / + a" erf 



a 



- a M \ l\/2crV n y I y^n - 1) 



where erf is as usually the error function and erfc its complementary with respect to 1. The 
second addendum in the brakets gives non vanishing contribuition only when x grows like n, 
therefore we can ignore it. What remains is a gaussian plus a correction that is significant when 
x is near to the perturbation. Infact defining 



erf a (x) = f dte-^^ 2 
Jo 



we have that 



erf CT (x) 



erf 



a 

so we can rewrite our correction as 



, x . i / s \/2vr|d| sign(x) ^ 2 j / larl In — 1 

(21) A* (a) = 1 1 Jt v ' e 1- erfc CT ' ! ! 



Case v = 2: In dimension two it is 



■ ^ir— 2n\B\ 2 (b~K,x) 
l a 2 (l — a) \x 



d a ~ 2?i \ = ril e 2n ( e 2n(n_1) — e~~^~ + o(l, x) ) , 



where o(l, x) is a function of x vanishing when n — > oo when x does not grow like n. Thus we 
can write 

. 2\B\\d\caa{d,x) {b-K, x) ( ^n-i 

(22) Agfa) = 1 " 1 , , v ' ' e ss — e M»-D - e ™ — + o(l, x) 

n|x| \ 

We notice that this term gives contribution only near the perturbation. Infact if x is fixed, i.e. 
not growing with n, A n (x) is a correction of the same order of diffusion. Moreover when x 

approaches to the origin it vanishes, since for large n e 2 ™(™-i) ~ 1 and 

^ re — 1 

lim j — j = 0. 

\x\ 

The correction disappears for larger x. If for instance x grows like yjri we easily see that A n (x) 
gives no contribution to the scale 1/n. 



Case v = 3: 
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The three dimensional case is qualitatively analogous to the two dimensional one. By an ele- 
mentaty change of variables e~ y = e~^ B x > x )/ 2na we can rewrite our integral as 

2n \ 3/2 r 2 . e~y y 2 



dy- 



2(„-l) (y_I__^J) 

, /2 (g^) V^FerfB (^) -2 v ^|x|e" i ^ 
= n C " ** 15-113/21x13 • 

up to corrections o (n 3 / 2 ) , and we notice that the limiting value of the last integral for x — > is 
|n 3 / 2 . In this way we can write at the leading order in n 

a3/ \ 1 COs(d,x)\d\ (g' 1 ^) / r— ( X \ / (.B~ l x,x) \ 

(23) Al(x) = ^ ^ ' e 2 n ^erf B ^J -2 v 0B rT |Ne 2 J. 

Finally we see that the correction term is relevant only when x is near the perturbation (though 
it vanishes at the perturbation point), but it gives higher order contributes in time to the asymp- 
totics when x is growing with n. 



These are expected results in high dimension. Of course a similar analysis can be pursued starting 
by (|18p for any value of the dimension v. 



Conclusions and outlooks 

In this paper we have studied a random walk on 1, u perturbed in one point. We have found 
corrections to the diffusion in the spirit of local limit theorem in the case of antisymmetric 
perturbation. 

Natural extension of the result is the symmetric case. As we mentioned, technical difficulties 
arise due to more involved form of the transition probability. Although we cannot show details 
about that, it seems clear that the behaviour of the correction should remain qualitatively the 
same. This is supported also by comparison to other known results [3], |S]. 

Another intersting perspective could be the analysis of the more general case of several points 
in which a perturbation is present, sparse, or all confined in a finite region. This would make a 
connection with the studies about obstacles in random walk in random environment f9|. 
We hope to report soon on these topics. 
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